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ABSTRACT 

In the region of Sgr B2 there are several condensations heated externally by nearby hot stars. 
Therefore H2O far-IR lines are expected to probe only an external low-density and high temperature 
section of these condensations, whereas millimeter-wave lines can penetrate deeper into them where 
the density is higher and T/j lower. We have conducted a study combining H2O lines in both spectral 
regions. First, Infrared Space Observatory observations of several H2O thermal lines seen in absorption 
toward Sgr B2(M) at a spectral resolution of ~35 km s _1 have been analyzed. Second, an IRAMSOm 
telescope map of the para-t^O 3i3 — 220 line at 183.31 GHz, seen in emission, has also been obtained 
and analyzed. The H2O lines seen in absorption are optically thick and are formed in the outermost 
gas of the condensations in front of the far-IR continuum sources. They probe a maximum visual 
extinction of ~5 to 10 mag. Radiative transfer models indicate that these lines are quite insensitive to 
temperature and gas density, and that IR photons from the dust play a dominant role in the excitation 
of the involved H2O rotational levels. In order to get the physical conditions of the absorbing gas 
we have also analyzed the CO emission toward Sgr B2(M). We conclude, based on the observed CO 
J=7-6 line at 806.65 GHz with the Caltech Submillimeter Observatory, and the lack of emission from 
the far-IR CO lines, that the gas density has to be lower than ~10 4 cm -3 . Using the values obtained 
for the kinetic temperature and gas density from OH, CO, and other molecular species, we derive a 
water column density of (9±3) x 10 16 cm" 2 in the absorbing gas. Hence, the water vapor abundance 
in this region, x(H20), is ~(l-2)xl0 -5 . The relatively low H2O/OH abundance ratio in the region, 
~2-4, is a signature of UV photon dominated surface layers traced by far-IR observations. As a 
consequence the temperature of the absorbing gas is high, Tk —300-500 K, which allows very efficient 
neutral-neutral reactions producing H2O and OH. On the other hand, the 183.31 GHz data provide a 
much better spatial and spectral resolution than the far-IR ISO data. This maser line allows to trace 
water deeper into the cloud, i.e., the inner, denser (n(iJ 2 )>10 5 ~ 6 cm" 3 ) and colder (Tfc~40 K) gas. 
The emission is very strong toward the cores. The derived water vapor abundance for this component 
is a fewxlO -7 . There is also moderate extended emission around Sgr B2 main condensations, a fact 
that supports the water vapor abundance derived from far-IR H2O lines for the outer gas. 
Subject headings: infrared: ISM: lines and bands — ISM: individual (Sgr B2) — ISM: molecules 



1. INTRODUCTION 

The determination of water abundance in space is 
a long standing problem in Astronomy. Theoretical 
models predict that water can be the most abundant 
species in warm m olecular clouds after H2 and CO 
ijNeufeld et al. 1995[) . Therefore, the determination of 
its spatial distribution and abundance contributes to a 
better knowledge of the chemical and physical processes 
that take place in the interstellar medium (ISM). 

Unfortunately, water is an abundant molecule in our 
atmosphere making particularly difficult the observation 
of its rotational lines and vibrational bands from Earth. 
Even so, some observations of water lines have been per- 

Electronic address: cerni@damir.iem.csic.es 

1 Based on observations with ISO, an ESA project with instru- 
ments funded by ESA Member States (especially the PI countries: 
France, Germany, the Netherlands and the United Kingdom) and 
with participation of ISAS and NASA. 

2 Present address: Laboratoire d'Etude du Rayonnement et de la 
Matiere, UMR 811, CNRS, Observatoire de Paris et Ecole Normale 
Superieure, 24 rue Lhomond, 75231 Paris Cedex 05, France 



formed from groun d-based and airborn e telescopes: the 
6i 6 -5 23 at 22 GHz llCheung et al. 1969(1 . the 5i 5 -4 22 at 
325 GHz llMenten et al. 19901) . the 10 TO -9^ at 321 GHz 
l|Menten et al. 1990br) andthe 3 i3 - 2 20 at 183.31 GHz 
(Watson et al. 1980; Cernicharo et al. 1990, 1994, 1996, 
1999a; Gonzalez- Alfonso et al., 1998). Due to the maser 
nature of these lines, the analysis and interpretation 
of the spectra is not obvious. Among these lines, only 
two have been used to map the extended emission of 
water vapo r in Orion: the 3ia — 2 ?n -hereafter the 183 
GHz line- (iCernich aro et a l. 19941) a nd the 5is — 422 
at 325 GHz (Cernic haro et al. 1999j) . However, al- 
though we know from ISO and SWAS observations that 
water is extended (~25'x25') in Sgr B2 (Cernicharo 
et al. 1997, Neufeld et al. 2003, Goicoechea et al. 
2004), little is known about its excitation conditions 
and its detailed spatial distribution. An alternative 
to indirectly estimate the water abundance in the 
Galact ic Center (GC) i s to use related spec ies such as 
HDO ITacn et al. 19901 IComito et al. 20031) or H,Q+ 
( Phill ips et al. 19921: iGoicoechea fc Cernicharo 2 001). 
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In none of these cases, the determination of x(H20) is 
straightforward. 

The ISO mission (|Kessler et al. 19961) . and spe- 
cially, the Long Wavelength Spectrometer, LWS, 
l)Clegg et al. 1996|). and the Short W avelength Spectrom- 
eter. SWS. l)de Graauw et al. 19 96). have provided a 
unique opportunity to observe several H2O lines in a 
great variety of astronomical environments. Never- 
theless, the majority of these observations were per- 
formed at the low spectral resolution of the grating mode 
(^1000 km s _1 ), which produces a critically strong di- 
lution in the search for molecular features in most ISM 
sources. Nevertheless, the Sgr B2 cloud has been ana- 
lyzed and studied in detail with the LWS / Fabry-Perot, 
which provided a velocity resolution of ^35 km s _1 
(Goicoechea et al. 2004, hereafter G04). 

Opposite to what is f ound toward othe r star forming 
regions such as Orion l|Cernicharo 1 999) . the observa- 
tions of the 2i2 — loi line at ~179.5 /an in Sgr B2 show 
that t he line appears in abso rption rather than in emis- 
sion ((Cernich aro et al. 1997^1 . Afterward, the launch 
of the Submillimeter Wave Astro nomy Satellite, SW AS, 
l|Melnick et al. 20001) . and ODIN l|Nordh et al. 2003J) al- 
lowed the observation of the lio — loi fundamental tran- 
sition of both H^O at 557 and YLfO at 548 GHz, first 
detected by the Kuive r Airborne Observatory, KAO, 
l)Zmuidzinas et al. 19951) . Although the velocity resolu- 
tion is ~1 km s _1 , the large beam of SWAS (~4') makes 
these observations more sensitive to the cold and less 
dense gas. SWAS observations have provided a reliable 
estimate of the water vapor abundance in the low exci- 
tation clouds located in th e line of sight toward Sgr B2 
(|Neufeld et al. 2000t 12001 . However, the fact that only 
the ground-state absorption line is detected makes dif- 
ficult a detailed study of water vapor excitation mecha- 
nisms in Sgr B2 itself. This is the most massive clou d 
in the Galaxy, with ~10 7 M Q , l|Lis fc Goldsmith 1990j) . 
and a paradigmatic object in the GC region as its ge- 
ometrical properties, physical conditions and chemical 
characteristics make it a miniature galactic nucleus with 
~15' extent (G04 and references therein). 

The main star-forming regions in Sgr B2 are located 
within three dust condensations, labeled as (N), (M) 
and (S). These condensations are embedded in a ~10 pc 
moderate-density (n(H 2 ) = 10 5 — 10 6 cm" 3 ) clou d 
l)Lis k Goldsmith 199ll: IHuttemeister et al. 1993|) . 
In addition, these structures are surrounded by 
lower density components of warm (Tfc>100 K) 
gas, hereafter the Sgr B2 warm envelope. These 
conditions have been derived mainly from ab- 
sorption observati on s of NH3 metastable lines 
(|Wilson et al. 19821 IHuttemeister et al. 1995]) . OH 
lines (Goicoechea fc Cernicharo, 2 001; 2002), and H 2 CO 
lines ( Ma rtin-Pintado et al. 1990j) in the radio domain. 
Nevertheless, the warm and low density gas is poorly 
traced by radio observations of other molecular species 
(generally excited by collisions in the denser regions and 
thus observed in emission). However, the warm envelope 
represents the strongest contribution to the absorption 
features produced b y man y light hydrides in the far-IR 
spectrum of Sgr B2 llcol . The origins of the observed 
rich chemistry and the heating mechanisms in the 
Sgr B2 warm envelope are a subject of intense debate. 
The matter is complicated due to the different observa- 



tional signatures to be integrated in the same picture: 
high temperatures derive d from NH3 absorption lines 
ijHuttemeister et al. 19951 iCeccarelli et al. 2002T) , fine 
structure emissio n from the photo-ionized and photo 
dissociated gas IIG04D. SiO and X-ray distribution 
l)Martin-Pintado et al. 20001) . etc. 

In all possible scenarios, water plays a significant 
role. Several mechanisms allow its formation and sur- 
vival in the warm envelope. Dissociative recombina- 
tion of H30 + leads to the production of H2O and OH. 
This processes depend on the specific Jh 2 o branch- 
ing ratio for the H 2 formation channel. Unfortu- 
nately, the determination of Jh 2 o with different ex- 
perimental proced ures has also yielded different values, 
from .^0=0.05 ( Williams et al. 1996]) to / ff2 o=0.25 
l|see Jensen et al. 20 00). while most chemical models 
have used /# 2 o~0.35. In addition, water could also be 
produced in the gas phase by the endothermic reaction: 

OH + H 2 ^ H 2 + H (1) 

However, the gas temperature must exceed ~300 K to 
overcome the activation barrier l)Neufeld et al. 19951) . At 
these temperatures, the reaction: 

O + H 2 -> OH + H (2) 

also contributes to the formation of OH. Therefore, H2O 
and OH column densities can be used to determine 
the role of the neutral-neutral reactions in their forma- 
tion/destruction routes. Still, the exact H 2 0/OH ratio 
will be determined by fn 2 o-, the temperature and by pho- 
todissociation processes if UV ra diation is present. As 
an example, Neufeld et al. l)2002|) studied a diffuse cloud 
(Go^l, n(H)=100 cm -3 ) toward W51 and showed that 
the presence of a warm gas component (Tfc>400 K) could 
explain the observed variations of the H 2 0/OH ratio re- 
spect to other diffuse clouds. 

Finally, high oxygen depletion onto water ice man- 
tles in dust grains could have taken place during the 
evolution of the cool gas in Sgr B2. Photodesorption 
and/or evaporation for dust temperatures above ~90 K, 
could release some water back into the gas phase enhanc- 
ing the H2 O abundance expe cted from pure gas-phase 
formation ijBergin et al. 2000T) . However, gas and dust 
are thermally decoupled in the outer layers of Sgr B2, 
where the d ust temperatures are sig nificantly lower, 
T d ~20-30 K llGordon et al. 1993HG04I). than gas tem- 
peratures, Tfc~300 K ijGoicoechea fc Cernicharo 20021 
ICeccarelli et al. 2002D . Thus, T^ seems too low to pro- 
duce significant evaporation of water ice mantles. There- 
fore, a detailed study of the far-IR H2O lines and of the 
183.31 GHz extended emission is needed to constrain the 
water abundance and the physical characteristics of the 
absorbing/emitting region. 

In this work we present and analyze the far-IR obser- 
vations of several thermal lines of water vapor toward 
Sgr B2(M) and the first map of the 183.31 GHz maser 
emission of para-H20 around Sgr B2 main condensa- 
tions. The layout of the paper is as follows: In Section[21 
we summarize the far-IR, submm and mm observations 
and data reduction. The spectra and maps are presented 
in Seed! Section 01 is devoted to the analysis of CO ob- 
servations (|4.1|l and water vapor observations Ij4.2|l with 
different radiative transfer methods. The main implica- 
tions of our work are discussed in Sec. [SJ where photo- 
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chemistry models for H2O and OH are also presented. A 
summary is given in Sec. El 

2. OBSERVATIONS AND DATA REDUCTION 

2.1. Far-IR observations 

Most pure rotational lines of water vapor that play a 
role in the radiative heating and cooling of the dense 
ISM appear in the Terahertz domain. Thus, airborne 
or satellite observations are needed to avoid the Earth's 
atmosphere blocking. In particular, many pure rota- 
tional lines of H2O a ppear in the far-I R coverage of 
the LWS spectromet er l|Clegg et al. 1 996) on board ISO 
ijKessler et al. 1996J) . We have used the LWS Fabry- 
Perot (FP) instrument to search for H2O and H^O lines 
toward Sgr B2(M). The LWS/FP spectral resolution is 
A/AA ~ 7000-1000 and has a circular aperture of about 
80" in diameter. The majority of detected lines have 
been observed in the time awarded to our ISO propos- 
als. However, an extensive inspection of the public ISO 
data base 3 has been carried out in order to examine and 
average all available water lines. Most of the water lines 
were present in the Astronomical Observation Template 
(AOT) LWS04 observations, which give a large spec- 
tral sampling, wavelength precision and S/N ratios. The 
number of scans in this mode was >12 depending on the 
expected absorption produced by the different species. 

These LWS products have been processed and com- 
pared through Off-Line-Processing (OLP) pipeline 6.0 
and 10.1 versions. There are no major differences ex- 
cept that recent pipelines produce <10% less absorption 
in some lines due to continuum level differences from 
one OLP to another. The data were analyzed using 
the ISO spectrometers data reduction package ISAP 4 . 
The mean FP continuum flux of each line deviates by 
<20% and this can be taken as the flux calibration error 
l|Swinvard et al. 19981 iGrv et al. 20f)l . After checking 
the continuum level of the AOT LWS01 observations, 
a polynomial baseline was fitted to each spectra and 
adopted as FP continuum level. 

2.2. Millimeter and submillimeter observations 

The water line at 183.31 GHz can be observed 
from high mountain-top sites under very dry conditions 
ijCernicharo et al~ 990). The observations in the direc- 
tion of Sgr B2 were achieved with the IRAM-30m tele- 
scope, with a half-power-beam-width (HPBW) of ~13", 
in December 1999 and January 2004. The source tran- 
sits at 24 degrees elevation at Pico Veleta, but even so 
the average atmospheric opacity over a bandwidth of 0.5 
GHz centered at 183.3 GHz was 1.7, i.e. an atmospheric 
water column of 0.9 mm along the line of sight. We 
used and SIS receiver designed to cover the band 130- 
184 GHz. In its higher end, the receiver temperature 
is about 160 K and the sideband rejection is at least 
18 dB. The backend used was a 512 two-pole filter with 
half-power widths and spacings equal to 1.0 MHz. With 
these conditions, and in this configuration, system tem- 
peratures were typically around 3000 K. The pointing 

3 The ISO data base observations (TDTs) presented in these 
work are 32201428, 32201429, 46201118, 46201123, 46900332, 
47600907, 47600908, 47600909, 47601001 and 49800301. 

4 ISAP is a joint development by the LWS and SWS Instrument 
Teams and Data Centers. Contributing institutes are CESR, IAS, 
IPAC, MPE, RAL and SRON. 



was check ed using the already k nown strong emission 
of W49N l)Cernicharo et al. 19901) . Once it was verified 
that the Sgr B2 emission was quite compact, the map was 
carried out in wobbler-switching mode in order to obtain 
very flat baselines. Observations were also tried in Jan- 
uary 2001, 2002 and 2004. In the 2001 run the weather 
was poor while in the last run it was possible to have a 
30 minutes window of observing time with good atmo- 
spheric transmission. These runs allowed us to check that 
the line profile in the position peak of Sgr B2(M) was the 
same of December 1999 and to confirm that the extent of 
the water emission toward this condensation was larger 
than the beam size. In the January 2004 observing run 
we also had very good atmospheric transmission to per- 
form observations around Sgr B2(N). However, the mea- 
sured line fluxes were lower than in previous observations 
by a constant factor. As the line profiles did not change, 
we scaled the 1999 observations to those the reference 
taken in the (0",+5") position in 2004 (see Fig.EJ. 

The CO 7-6 (806 GHz) observations were performed 
on March 10, 2002 with the 10.4 m Caltech Submillime- 
ter Observatory (CSO) located at the summit of Mauna 
Kea (Hawaii). The receiver is a helium-cooled SIS mixer 
operating in double-sideband mode (DSB) providing an 
instantaneous bandwidth of 0.95 GHz and designed to 
fully cover the 780-910 GHz atmospheric window (see 
Kooi et al., 2000). The pointing was checked using a 
strong CO point-like source in the nearby W28 molecular 
cloud, and was kept within 3" accuracy for a HPBW of 
~10". Two different acousto-optic spectrometer (AOS) 
backends were used with 2048 and 1024 channels respec- 
tively for a total spectral coverage of 1.5 and 0.5 GHz 
in each case. The zenith atmospheric water vapor col- 
umn was ~0.5 mm during the observations, resulting on 
system temperatures ranging from 6000 to 10000 K de- 
pending on the receiver and the elevation. CO emission 
is very extended around Sgr B2, so we had to perform 
position switched scans setting the off position 1 degree 
away in azimuth. An additional 5% increase had to be 
applied to T^ in order to second-order correct the stan- 
dard chopper-wheel" calibration method at high frequen- 
cies with large atmospheric opacities (Pardo et al. 2005). 

3. RESULTS 

The far-IR spectrum of Sgr B2(M) is dominated by the 
absor ption produced by NH3, OH and H2O rotational 
lines (G04). The detected far-IR lines of water vapor 
are shown in Fig. ^ H2O spectroscopical and observa- 
tional data are tabulated in Tabled Except the 2i2 — lot 
ground state line of 0-H2O at ^179.5 /jm (~1669.9 GHz), 
all water lines have a similar profile and are centered at 
Sgr B2(M) velocities. The ~179.5 /xm line is saturated 
and absorbs from -150 to +100 km s _1 , which there- 
fore includes the water vapor located in the foreground 
gas toward the GC and the warm gas around Sgr B2(M). 
The widespread absorption produced by the 2i2 — lni line 
has been previously presented in Ccrnicharo et al. (1997) 
and enlarged in G04 (9'x27'), while the lio — lot absorp- 
tion has been mapped (26'xl9') by SWAS (Neufeld et 
al. 2003). These observations probe that low excitation 
H2O is present in the clouds intersecting the line of sight 
toward large areas of Sgr B2. 

The average velocity of all H2O lines observed with the 
ISO/LWS-FP is +60±5 km s" 1 , in agreement with the 
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Fig. 1. — ISO/LWS/FP observations of H3, 6 and HJ, s O toward Sgr B2(M) and rotational energy diagram of ortho-H3, 6 0. The ordinate 
scale corresponds to F/F c and the abscissa to the wavelength in microns. All the far-IR water lines appear in absorption. The emission line 
corresponds to the 3i3— 220 maser transition of para— H^O at 183.31 GHz (observed with the IRAM-30m telescope at Pico Veleta, Spain). 



velocity of other related oxygen species such as H30 + 
or OH (Goicoechea & Cernicharo, 2001; 2002). Tak- 
ing into account the wavelength calibration error of the 
LWS/FP instrument, this velocity is compatible with the 
expected +65 km s" 1 clo ud seen in radio observations 
ijHuttemeister et al. 19951) . We note, however, that ve- 



locities close to +60 km s -1 are also associated with gas 
surrounding most Sgr B2 continuum sources, so that the 
bulk of the H2O absorption can arise from them. 

Possible overlapping with other molecular species oc- 
curs at some wavelengths. In particular, the 0-H2O 
4 3 2-4 2 3 line is blended with HF J=2-l at 121.697 /im 
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Fig. 2.— Map of the 3i 3 -2 2 o line of para-H^ 6 at 183.31 GHz around Sgr B2 main condensations. Contours in K km s arc indicated 
in the figure. Different line profiles at significant positions of the map are also shown. The intensity scale is in and the abscissa is the 
LSR velocity in km s — x . The map is centered at the position of Sgr B2(M). The spectra were obtained in two different runs in Dec. 1999 
and Jan. 2004. The fluxes were scaled using the (0" ,5" ) position as reference. Note that the line profiles in that reference position did not 
change in more than 4 years. 



(jNeufeld et al. 1997ft and the o-H^ 8 2i 2 -l i line at 
181.053 /xm has a small contribu tion from H3O 4 " Q(l,l) 
ijGoicoechea fc Cernicharo 2001ft . 

Observations of the 183.31 GHz H 2 line (E„~205 K) 
toward Sgr B2 are presented in Fig. The emission ap- 
pears at the LSR velocities of Sgr B2 with no contribu- 
tion from the line of sight clouds. The emission appears 
at least in the 50-75 km s _1 range. Note the different line 
shapes and intensities of the 183.31 GHz emission for po- 
sitions in front and around the main condensations (here 
the line appears much wider). 

Finally, Fig. shows the CO J=7-6 line (E„~150 K) 
observed toward Sgr B2(M). This CO line shows a strong 
self-absorption at sim70 km s _1 (the velocity of the 
183.31 GHz H2O line peak at his position) so that the 



profile peaks at ~ +55 and +85 km s . Similar patterns 
are shown by lower- J CO lines. 

4. ANALYSIS AND DISCUSSION 

4.1. Carbon Monoxide 

The warm gas present in the outer layers of Sgr B2(M) 
might be expected to radiate in high-J CO lines. The 
CO J=14— 13 transition at 185.999 /im is the one with 
the lowest energy level (E u ~581 K) within the range 
of ISO/LWS detectors. However, we have not de- 
tected any emission/absorption from CO (3<r limits are 
<2x 10~ 18 W cm" 2 ) in the ISO/LWS spectra toward the 
Sgr B2 region |G0j), with both grating and FP. CO spec- 
troscopical data and line flux upper limits are tabulated 
in Table |3 Nevertheless, recent studies of the large scale 
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Fig. 3.— CSO observation of the CO J=7-6 line at 806 GHz 
toward Sgr B2(M). 



CO J— 7-6 emission toward the GC with the Antarc- 
tic Submillimeter Telescope and Remote Observatory 
(AST/RO) have shown that the emissi on is concentrate d 
toward Sgr B and Sgr A complexes ijKim et al. 2 002'). 
Fig. shows higher spectral/angular resolution observa- 
tions of this line taken with the CSO telescope toward the 
Sgr B2(M) position. Hence, it seems that at a given J up 
level, the rotational CO line emission disappears from 
the Sgr B2 spectrum. 

We have performed nonlocal radiative transfer calcula- 
tions to try to reproduce the lack of high-J CO lines in 
the far-IR spectrum of Sgr B2(M), and to help constrain- 
ing the physical parameters needed to model the water 
vapor absorption. In Fig. 0] we show the predictions of 
a nonlocal model for several high-J transitions of CO. 
The nonlocal model used here is an adaptation of the 
radiative trans fer co de developed by Gonzalez-Alfonso 
& Cerni charo ll!993f) with the inclusion of dust in the 
transfer ijCernicharo et al. 200 0"). We have implemented 
it for high-J levels of CO (see Sec. 14.2.21 for further 
details on the specific model developed for Sgr B2(M)). 
The collisional rates have been taken from Flower (2001). 
The dust continuum emission has been modeled follow- 
ing Goicoechea and Cernicharo (2002) and G04 (see sec- 
tion 4.2.1). We have run an array of models for CO 
abundances 5 of 10 -4 and 10 -5 with varied from 100 
to 500 K and n(H 2 ) from 5.0xl0 3 to 6.4xl0 5 cm" 3 . 
From these results it is clear that in order to suppress the 
far-IR CO emission, and to match the CO J=7-6 line 
emission, low H2 density is required. In particular, if the 
CO J— 7-6 line arises from a layer of gas at T^^lOO K 
the model with n(ff2)=2x 10 4 cm -3 correctly reproduces 
the absence of far-IR lines. However, if it arises from the 
same outer layer of warm OH (Tfc~300 K ) detected in 
the far IR (Goicoechea & Cernicharo 2002), the limit for 
the density will be <10 4 cm -3 . In any case, the lack of 
high-J CO lines is an evidence of the low density gas 
located in front of the far-IR continuum source. In the 
following section we will constrain the physical parame- 
ters of this layer by studying the available water vapor 
lines. 

5 An upper limit to the 13 CO fractional abundance in Sgr B2 
of 10 — 6 was found by Lis & Goldsmith (i989j). According to 
12 C 18 0/ 13 C 18 0~25 in Sgr B2 ILanger Ik Penzias 19901) . the 12 CO 



4.2. Water Vapor 

The main observational result of this work is that the 
far-IR water lines toward Sgr B2(M) appear in absorp- 
tion, while the 183.31 GHz line is seen in emission around 
and in the main condensations. The fact that the contin- 
uum emission in the far-IR is optically thick l|G04f) in- 
dicates that the H2O absorption lines arise from regions 
where the excitation temperatures (T ex ) are smaller than 
the dust temperatures inferred from the continuum emis- 
sion. However, the problem of the H2O line excitation 
toward Sgr B2 is not straightforward. Apart from self- 
absorption and the possible excitation by collisions with 
molecules (e.g. H2), atoms (e.g. He) and e~ (if the ion- 
ization fraction is significant), the level population can 
be primarily determined by the thermal emission of dust, 
the role of which is essential in the excitation of molecules 
such as H2O or OH which have many rotational lines in 
the far- and mid-IR. This represents a major difference 
with respect to the excitation treatment of molecules ob- 
served in the radio domain where, generally, one can ne- 
glect the excitation by dust photons. If the excitation 
is dominated by IR photons, the two transitions arising 
from the ground levels of ortho-water, 2i2-loi and lio - 
loi, will determine how the higher energy levels will be 
populated. Even if collisions are important, the pres- 
ence of an optically thick far-IR continuum will strongly 
affect the T ex of water and, thus, it must be carefully 
taken into account in the models. In the case of Sgr B2, 
this means that the external dust layers of the cloud will 
absorb the possible water line emission from the inner 
regions. Knowledge of the geometry of the region to be 
modeled and the relative filling factors of the dust and 
gas in the beam of the LWS instrument, are also im- 
portant for the models. For this reason, high angular 
resolution ground-based observations of the 183.31 GHz 
water line are particularly important to model the water 
vapor radiative transfer in Sgr B2. 

4.2.1. Large Velocity Gradient modeling 

In this section we analyze the H2O observations with 
a multi-molecule Large Velocity Gradient (LVG) model. 
For Sgr B2(M) we have adopted a spherical geometry 
with two components (see Fig.0: a uniform continuum 
core with a diameter of ~23" (~1 pc for a distance of 
8.5 kpc) and a shell of variable thickness and distance to 
the core. The presence of an external shell of molecu- 
lar gas (not resolved by the ISO/LWS beam) surround- 
ing a central cond ensation is indicated by the an alysis 
of the far -IR OH ijGoic j^hj^ fc Cernicharo 200*21) and 
NH3 lines ijCeccarelli et al. 20 02) and it is also supported 
by H2O 183.31 GHz line observations (Fig.|5J|. In partic- 
ular Goicoechea & Cernicharo (2002) found an angular 
size of ~42" for the OH envelope. The bulk of the H 2 
absorption lines must arise from the OH layers, or inside 
them. Following the OH model geometry and Figs.|5]and 
03 a total size of ~30" (~1.2 pc) for the core+shell cloud 
has been adopted. The core is considered as a gray-body 
with an opacity at 80 (im of 2.5, with a dust opacity law 
given by t\ = T8o[80/A(/xm)], and a dust temperature 
of 30 K. These values are consistent with the color tem- 
peratures and dust emissivities derived from the analy- 
sis of the ISO/LWS continuum observations at the same 

abundance would be <2.5xl0 — 5 . 



Water Vapor around Sgr B2 



7 




10 15 20 25 

J, 

Fig. 4. — Nonlocal models of CO rotational lines toward Sgr B2. Velocity-integrated line fluxes (in 10 — 17 W cm -2 ) are shown as a 
function of J U p, the CO upper level rotational number. From J=14-13, CO lines are accessible to the ISO/LWS (vertical dashed line). No 
CO line has been detected. 3<r upper limits for CO line fluxes are shown (filled squares). Two different abundances have been considered, 
x(CO)=10 — 4 (upper panels) and x(CO)=10 -5 (lower panels). Models with different kinetic temperatures (100 to 500 K) are shown in each 
panel for different H 2 volume densities: 5.0xl0 3 , l.OxlO 4 , 2.0xl0 4 , 4.0xl0 4 , S.OxlO 4 , 1.6xl0 5 , 3.2x10 s , 3.2xl0 5 and 6.4x10 s cm" 3 . 



wavelengths of the detected far-IR H 2 lines ljG04f) . Al- 
though dust temperatures can be slightly larger (T<j~60- 
80 K) inside the cloud (from the analysis of the millimeter 
continuum emission) or slightly lower (Td~20 K) in the 
outer and colder diffuse layers (from the analysis of ex- 
tended IRAS continuum emission; Gordon et al. 1993), 
we judged T^— 30 K as the most representative value for 
the dust grains coexisting with the molecular species de- 
tected in the far-IR. The bulk of the continuum emission 
detected by ISO can be hardly fitted with dust temper- 
atures below 30 K. The emission of higher temperature 
dust arising from the innermost regions of Sgr B2 is hid- 
den in the FIR by the huge amount of foregroung gas 
and colder dust. 

For these models we have considered an ortho-B^O 
column density of 1.8xl0 16 cm -2 and a para-B^O col- 
umn density of 0.6xl0 16 cm~ 2 . These are lower lim- 
its suggested by the nonlocal radiative transfer mod- 
els (see below). The LVG model computes the sta- 
tistical equilibrium population of the rotational levels 
for ortho-B20 and para-B20 independently. The colli- 
sional rates were sc aled from those of B2 O-Be collisions 
HGreeu et al. 199.31). 

Several LVG computations for some selected ortho- 
and para-B20 transitions in the TBz domain are shown 
in Fig. |SJ The excitation temperature T ex of each tran- 
sition is shown in each panel as a function of T^ and 
n(i?2). The thick contour corresponds to the equivalent 
temperature of the continuum core. Therefore, to the left 
of this contour, water lines appear in absorption. From 
the LVG models it is clear that high density and tem- 
perature are required to observe the far-IR water lines 
in emission. For the range of densities implied by the 
CO observations (<10 4 cm -3 ), the B2O lines observed 
by the ISO/LWS are correctly predicted in absorption. 



The temperature of the absorbing layer is, however, more 
difficult to estimate because the possible solutions for a 
given (low) density model are not very sensitive to tem- 
perature variations, as indicated by the smooth change 
of T ex for constant density as the kinetic temperature 
goes from 500 to 100 K. LVG models predict that, in 
addition to dust photons, collisions play a role in the ex- 
citation of the lowest B2O rotational levels. The higher 
energy levels are pumped from the lowest levels by ab- 
sorption of far-IR photons. In the case of the ortho-H^O 
lio-loi line observed by SWAS (Neufeld et al. 2000), 
only moderate densities (>5xl0 4 cm -3 ) are needed to 
observe the line in emission. This is the case of the ex- 
tended emission in the lio-loi line observed in the 180 pc 
Molecular Ring around the GC between i>ls.r^+80 and 
120 km s -1 (Neufeld et al. 2003). For the density con- 
ditions derived for the Sgr B2 velocity range and for the 
line of sight clouds, the 557 GBz line is correctly pre- 
dicted in absorption. A similar behavior is expected for 
the para-B20 ln-loo line at ~269.3 /um (^1113 GBz) 
that will be observed by future heterodyne instruments 
such as BIFI/Berschel. 

Among all the moderate excitation para-B20 lines 
(E;<450 K), only the 3i3-2 2 o appears in the mm do- 
main. Contrary to other B2O maser lines accessible from 
ground-based telescopes, relatively low T^ and density 
are required for the 3i3-22o line inversion. These condi- 
tions allowed the first detection of extended water emis- 
sion in Orion (Cernicharo et al. 1994), and a strong de- 
pendence of the emission with T& was revealed. Bence, 
the 183.31 GBz line could be an excellent tracer of the 
warm gas in molecular clouds. 

We have used the same LVG model to analyze the 
183.31 GBz line inversion mechanism that produces ex- 
tended emission in Sgr B2 (Fig.[5J). Fig.[|)]shows different 
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Fig. 5. — Large Velocity Gradient excitation models of different kinetic temperature and density for several water vapor rotational lines. 
The column density of ortho-H20 is 1.8xl0 16 cm~ 2 and that of para-H20 is 0.6xl0 16 cm -2 . In each panel, the H2O transition and 
wavelength/frequency is labeled. The thick contour corresponds to the equivalent temperature of the continuum source (T c ). Hence, to 
the left of this contour, lines will be in absorption (T e:E <T c ). 



excitation models for the 3i3-2 2 o transition. Since a frac- 
tion of the 183.31 GHz line emission may arise from the 
warm envelope in front of the far-IR continuum source, 
models with and without a central continuum source (the 
same described at the beginning of the section) have been 
considered. For a given T& (from 20 to 500 K), each panel 
shows computations of different para-H 2 column den- 
sities as a function of Tb (in K) and H2 volume density. 
For low temperatures (Tfc<40 K) and moderate densi- 
ties (<10 5 cm -3 ), the line will be observable for large 
.ZV(p-H 2 0) values only if a continuum source is present. 
Due to the minor role played by collisions in the pump- 



ing of the rotational levels at these temperatures, the line 
intensity is almost independent of the density. Neverthe- 
less, as the temperature increases, collisions start to be 
significant, and the expected intensity of the 183.31 GHz 
line becomes less sensitive to the models with or without 
the continuum source. 

The 183.31 GHz emission is produced by the popula- 
tion inversion of the 3i3 and 2 2 o levels. This is due to 
the different rates (far from the thcrmalization densities) 
at which both levels can be populated. The 3i3 level is 
radiatively connected with the 2 20 and 2 02 lower energy 
levels. The A^- coefficient for the 3i3~2 02 at ^138.5 /tm 
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Fig. 6. — Large Velocity Gradient excitation models for the 3i3-220 maser transition of para-Hj O at 183.31 GHz. Models with 
and without a central continuum source (described as a gray-body with an opacity of 2.5 at 80 /jm, with a dust opacity law given by 
Tx = T8o[80/A(/xm)], and a dust temperature of 30 K) have been considered. For low temperature and low densities, infrared pumping seems 
to be efficient to increase the emerging flux of the 183.3 GHz line. In addition the maser amplifies the continuum. Hence, the resulting flux 
depends on the assumed opacity at 80 )im and on the opacity law. For higher kinetic temperatures and densities collisional pumping of 
the maser dominates and the difference between the two cases, with and without central continuum source, is less important. For a given 
kinetic temperature (20 to 500 K), the expected Tg as a function of the H2 density is shown in each panel for different column densities 
of para— H^O. The dashed line indicates the observed brightness temperature in the extended emission. Toward the main condensations, 
brightness temperatures of ~100 K are reached. 
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Fig. 7. — Sketch of the water model. The dimensions of the OH 
envelope are also shown (Goicoechea & Cernicharo 2002). 



is — 10 4 times larger than the Ay of the 3i3~22o transi- 
tion. The ~138.5 /im line is seen in strong absorption 
in the ISO/LWS spectra (see Fig. HJ. The 2 20 level is 
radiatively connected with the 2n and In lower energy 
levels. The 220-ln at —100.9 /im produces the strongest 
absorption of para-EbO in the far-IR (see Fig. QJ. On 
the other hand, the 3i3 level is radiatively connected with 
the 322, 4 2 2 and 4 4 higher energy levels, while the 2 2 o 
level is connected with the 331 (see the —67.1 /xm line 
in Fig. Q and with the 3i3 higher energy levels. The 
different radiative pathways and rates (for n«n cr ) at 
which the 3i3 and 2 2 o levels can be populated produce 
the inversion. LVG models (see Fig. [SJ show that this 
mechanism can be efficient to produce 183.31 GHz line 
emission in regions of relatively low density. 

The observed extended emission at 183.31 GHz has a 
brightness temperature of 10 K in average (see Fig. |5J. 
Assuming that it arises from the low density regions in 
which the the far-IR continuum sources are embedded, 
with T fc =300-500 K and n(iJ 2 )-10 4 cm" 3 , we derive a 
value for A^(p-H 2 0) of >5xl0 17 cm~ 2 . This value is 
higher than the one derived from the far-IR water lines 
(see sect. 14. 2.211 possibly because the 183.31 GHz line 
penetrates deeper into this dusty environment. Toward 
the main condensation, the bulk of the emission seems 
to arise from the cold (T/. —40 K) and dense gas. Under 
these conditions, a brightness temperature of the line of 
—100 K would translate into a P-H2O column density of 
— 10 19 cm -2 . However, an important contribution could 
come from the embedded high temperature and high den- 
sity core condensations where the column density should 
also be much larger. 

4.2.2. Non-Local radiative transfer models 

To take into account the radiative coupling between 
regions with different physical and/or excitation condi- 
tions, the radiative transfer has to be treated with nonlo- 
cal techniques, more sophisticated than the LVG approx- 
imation. We have adapted the radiative transfer used in 
Sec. 5.1 for CO, to the ground vibrational states of ortho- 
and para-H20 respectively. The model includes all the 
water rotational levels with transitions between 40 /im 
and 183.31 GHz. The level population is computed in 



statistical equilibrium considering collisional excitation 
and radiative excitation by line and continuum photons. 
This is computed consistently assuming that the water 
molecules and the dust grains are coexistent. The ge- 
ometry, core+shell dimensions, are the same of that con- 
sidered in the LVG models. The shell was divided in 41 
layers. The central dust condensation has been modeled 
with identical parameters than in section 4.2.1. We have 
considered 14 rotational levels (E„ < 608 K) of ortho- 
water for model with T& below 100 K , and up to 30 
rotational levels (E up < 1290 K) for the higher tempera- 
ture models. We have adopted a turbulence velocity of 8 
km s _1 . The continuum radiation fiel d has been treated 
as in Gonzalez- Alfonso & Cernicharo l)1993f) by consider- 
ing a spectral range of 70 km s _1 centered on each water 
transition. Collisional rates have been taken from Green 
et al. (1993). Obviously all the water transitions in the 
continuum core are thermalized to the dust temperature 
due to the large opacity in the far-IR. The computed line 
profiles are a result of convolving the brightness temper- 
ature with the angular resolution of the LWS detectors 
(~80"). The resulting spectral resolution of the syn- 
thetic water lines is 1 km s _1 , as no convolution with 
the spectral resolution has been performed. To test the 
sensitivity of the model to the physical parameters, mod- 
els were computed with A(H 2 0) and Tj, of 1.8xl0 16 , 
9xl0 16 and 4.5 xlO 17 cm ~ 2 , and 40, 100, 200, 300 and 
500 K respectively, while densities were increased from 
5xl0 3 to 6.4xl0 5 cm ~ 3 multiplying by 2 in each step. 
The different nonlocal radiative transfer models for the 
first two H2O column densities are shown in Figs. |S] and 
Models for other column densities have been also ran. 
However, the observed absorption depth is not repro- 
duced for column densities below 10 16 cm ~ 2 . In the 
A(H2O)=1.8xl0 17 cm ~ 2 case too many lines would be 
in emission, contrary to the observations. 

The main problem to interpret the H2O absorption to- 
ward Sgr B2 arises from the large opacities of the far-IR 
water lines, — 10 3 , and even — 10 4 for the ortho-H 2 2i 2 - 
loi line at —179.5 /im. Under these conditions, many 
weak water lines have to be detected to constrain the 
physical conditions and the column density. Therefore, 
we strength that it is difficult to obtain physical param- 
eters from the observation of few far-IR water rotational 
lines. In addition, as the radiative excitation dominates 
the population of the far-IR levels, lines with a weak de- 
pendence on the dust excitation should be investigated. 
Taking into account the 14 water lines detected in the 
far-IR, the nonlocal results imply that models are not 
very sensitive to the temperature and that the only in- 
dication about the column density has to be searched in 
weak far-IR H2O lines or in the 183.31 GHz line (see 
Sec. 6.1.2). Even so, the far-IR absorption arises in the 
low density external layers of gas, while it is very likely 
that the 183.31 GHz line could have an important contri- 
bution from inner and denser regions. Another complica- 
tion arises from the fact that for models of low A(H 2 0) 
values and Tj,<200 K, it is also difficult to distinguish 
between different H 2 densities (see Fig. |S). Limits to 
Tfc and n(H2) have to be searched in weak lines below 
70 /im. 

The models for high water column densities (Fig. 
predict absorption lines in the LWS range at —56.3, 
-57.6, -58.7, -78.7, -99.5, -125.4 and -136.5 /mi. At 
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Fig. 8. — Results from selected nonlocal models of H2O rotational lines in Sgr B2(M). Results for different gas temperatures (from 40 to 
500 K) and H 2 densities (from 5.0xl0 3 to 1.6xl0 5 cm -3 ) are shown. Note that for each model, the water abundance changes accordingly 
with n(H2) in order to keep constant the water vapor column density (in this case, 1.8xl0 16 cm -2 ). For the first nine lines, the frequency 
is given in GHz; for the rest, the wavelength in fim is indicated. The velocity scale is shown at the bottom of the 113.5 /im line. 



the spectral resolution and sensitivity of the LWS/FP, 
none of these lines has been detected. This implies that 
iV(H2O)<4.5T0 17 cm -2 toward the warm envelope of 
Sgr B2. The ~136.5 //m line (also predicted by the mod- 
els with large column densities) is contaminated by the 



absorption produced by the C3 Q(S) rovibrational line, 
whic h is also predicted by th e models of tri-atomic car- 
bon ijCernicharo et al. 200Cf) . Some of these lines are 
predicted (even in emission) by the models with large 
iV(H 2 0). Models with iV(H 2 O)=1.8xl0 16 cm" 2 are con- 
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Fig. 9. — As in Fig, ISlbut for a water vapor column density of 9.0xl0 16 cm 
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sistent with ISO detections and upper limits to other 
ISO lines. Only the ~67.1 /mi line is weaker than the 
observations. Therefore, models shown in Fig. [S] give 
a lower limit to the water vapor column density in the 
outer and warm (300-500 K) envelope. Taking into ac- 
count the difficulties implied by the H2O modeling in 
the far-IR, we found that (9 ± 3) x 10 16 cm -2 is the best 
H 2 column density to fit the ISO observations (see Fig- 
ure H5J. The CO analysis (Sec. 5), and the studi es in 
the far-IR of OH iGoicoechea fc Cernicharo 2002(1 and 
the ammonia lines l)Ce^clu"eTIT"et"lt!"" 20021) and our CO 
data, also support that the water absorption lines arise 
in the warm and low density (<10 4 cm" 3 ) layer in front 
of Sgr B2(M). The H2O column density derived from 
ISO observations is below the lower limit of A(H20) 
estimated for the 183.31 GHz line for this component 
(see Sec. 14. 2. If) as it is likely that an important frac- 
tion of the 183.31 GHz emission arises from the inner 
and denser regions closer to Sgr B2 main cores or even 
from them, as pointed out above. This component, with 
a mean H2 density of >10 5 " 6 cm" 3 , has been traced 
by the NH3 non-metasta ble emission lines, T,t~100 K, 
( Huttcmci ster et al. 1993ft and by the emission of HC3N 
rotational lines, T fe ~20-40 K, ijLis fc Goldsmith 199ljl . 
The determination of the temperature from mm emis- 
sion lines is also complicated in these regions completely 
obscured to ISO observations. Therefore, the observed 
differences in the 3 13-202 line-intensity and -shape can 
be a combination of A(H20) and/or T& variations across 
the region. 

The results presented in this section have been com- 
pared to those obtained from another radiative transfer 
based on a different approach (Asensio Ramos & Tru- 
jillo Bueno 2003; Asensio Ramos & Trujillo Bueno 2006). 
This code is a generalization to spherical geometry of the 
very fast iterative Multilevel Gauss-Seidel (MUGA) and 
Multilevel Successive Overrelaxation (MUSOR) methods 
developed by Trujillo Bueno & Fabiani Bendicho (1995) 
for the case of a two-level atom and generalized by Fabi- 
ani Bendicho, Trujillo Bueno & Auer (1997) to multilevel 
atoms in Cartesian geometries. The code also allows 
the application of the standard Multilevel Accelerated 
A-iteration (MALI) (Olson, Auer & Buchler 1986). The 
angular information for the calculation of the mean in- 
tensity is obtained by solving the radiative transfer equa- 
tion along its characteristic curves (straight trajectories 
with constant impact parameter) with the aid of the 
short-characteristics formal solver with parabolic preci- 
sion (Kunasz & Auer 1988) . The statistical equilibrium 
equations are linearized with the aid of the precondi- 
tioning scheme developed by Rybicki & Hummer (1991, 
1992) with the introduction of an approximate A oper- 
ator that can be efficiently obtained in the framework 
of the short-characteristics technique. The convergence 
rate of the MUGA and MUSOR schemes is equivalent 
to that obtained with the introduction of a nonlinear A 
operator, with the advantage of not being necessary nei- 
ther to calculate nor to invert such nonlinear operator. 
Interestingly, the time per iteration is similar to that ob- 
tained for the standard A-iteration or the MALI method. 
The computing time for the MUGA scheme is reduced 
by a factor 4 with respect to MALI, while the MUSOR 
scheme leads to an order of magnitude of improvement 
in the total computing time with respect to MALI. 



The calculations have been performed with the same 
geometry, the same molecular data (collisional and ra- 
diative transitions), and the same physical conditions of 
the previously described nonlocal code. The emerging 
line profiles from both codes are very similar, with dif- 
ferences below 2-3 %. Both codes predict lines in ab- 
sorption/emission for the same physical conditions, with 
identical spectral shapes in the cases where re-emission 
is found in the line wings. This test of consistency allows 
us to be very confident in the results presented in Figs. [HI 
and since they were obtained with numerical methods 
based on completely different approaches. 

5. THE WARM ENVELOPE: SHOCKS, PDRS OR XDRS ? 

The H2O (OH) column densities derived from ISO 
observations in the warm envelope are within an or- 
der of magnitude of the sh ock model predicti ons for 
the same low density region l)Flower et al. 1995[) . Nev- 
ertheless, C-shocks have been invoked to explain the 
heating of large amounts of war m gas in the GC 
and also in the warm en velope ijWilson et al. 198*21: 
iMartm-Pintado et al. 1 99^1 . In addition, the shock 
models of l|Floweret"'aL 1995|) correctly reproduce the 
observed column densities of A-bearing species such as 
NH2 and NH3 that c ould b e ultimately related with the 
dust grain chemistry l)O04f) . 

A qualitative explanation came from the observation 
of large OH column densities (H2 0/OH~2-4) in the 
warm envelope. According to Goicoechea & Cernicharo 
(2002), if a far-UV radiation field illuminates the 
outer regions of the cloud, water molecules could be 
photodissociated producing an enhancement of the OH 
abundance expected only from its formation through 
H3O" 1 " dissociative recombination. The presence of such 
an extended far-UV radiation field with 67o~10 3 -10 4 is 
inferred from the [Qui ], [N lll], [Nil], [C11], and [Ol] 
extended line emission (IG04T) . Therefore, an important 
contribution from H2O photodissociation may explain 
the observed H2O/OH abundance ratio in the external 
layers of the envelope. 

To investigate in more detail the O-chemistry in the 
envelope, we have run several photochemistry calcula- 
tions that include depth-dependent photodissociation, 
HsO" 1 " dissociative recombination and neutral-neutral re- 
actions. We have used the latest version of th e pub lic 
and available PDR model 6 of Le Bourlot et al. (1993). 
The model does not include oxygen grain surface chem- 
istry, but it is consistent with the low densities and 
hig h tem peratures found in Sgr B2 envelope. Follow- 
ing l)G04|h we assume a Go=5xl0 3 radiation field and a 
7i#=n(H)-t-2n(H2)=5xl0 3 cm" 3 density. Assuming that 
the outer gas layers of Sgr B2 envelope are directly illumi- 
nated by such a radiation field, the model solves the UV 
transfer and the chemistry up to Ay=20 mag. As a refer- 
ence model we solve the thermal balance explicitly for an 
initial gas temperature of 500 K. The resulting H2O and 
OH column densities and the H2O/OH ratio are shown 
in Fig. ^]| as a function of the visual extinction through 
the cloud. This pure PDR can only maintain the tem- 
perature above ^100 K in the first 2 magnitudes of the 
cloud, where the UV radiation efficiently photodissoci- 

6 http://aristote.obspm.fr/MIS/ 
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Fig. 10. — Selected PDR models in which Go=5xl0 3 and 
n^=5xl0 3 cm -3 . Upper panel: H2O and OH column densities 
as a function of the visual extinction for an isothermal model with 
T gas =500 K and for a model were the full thermal balance is solved 
at each depth. Lower panel: H2O/OH column density ratio as a 
function of the visual extinction for different isothermal models 
with T sas =500, 300, 200 and 50 K and for full thermal balance. 



ates H 2 to form OH. This results in a low H 2 0/OH 
ratio. Inside the cloud, the H2O and OH production is 
rapidly dominated by HsO + dissociative recombination 
and the H2O/OH ratio tends to a constant value that 
basically depends on the assumed branching ratio fH 2 o 
(we have taken /# 2 o=0.25). 

At least two different scenarios can produce large col- 
umn densities of gas at high TV The first one is the 
presence of several clumpy PDRs within ISO's beam. 
In this scenario, far-IR H2O and OH observations will 
only trace the PDR-clump surfaces. The second one is 
the presence of shocks. It is very likely that widespread 
low velocity shocks can locally heat the gas within the 
envelope and preserve a fraction o f the molecular ga s 
with temperatures up to ^500 K ijFlower et al. 19951. 
but the same is possible in the clumpy PDR scenario. 
With these temperatures neutral-neutral reactions play 
an important role in the chemistry. To simulate this 
situation, we have run several iso-thermal models with 
temperatures ranging from 500 to 50 K. Some selected 
results are also shown in Fig. ^3 From these models 
it is clear that both iV(H 2 0) and A^(OH) are clearly 
enhanced by neutral-neutral reactions (Fig. 1101 upper 
panel). If the temperature is high, these models re- 
produce the observational values much better. Hence, 
a significant fraction of warm gas seems to be needed to 
reproduce the large H 2 and OH column densities ob- 
served toward Sgr B2. This conclusion agrees with the 
large te mperatures derived from OH ob servations in the 
far-IR l)Goicoechea fc Cernicharo 20021) . 

The gas temperature determines the contribution of 
neutral-neutral reactions (1) and (2). As the temper- 
ature of the gas increases, reaction (1) contributes to a 
enhancement of the H 2 0/OH ratio fFig. HUI lower panel). 
Note that a moderate increase of the radiation field (to 



Go=10 4 ) shifts the H 2 0/OH curve to larger Av inside 
the cloud until photodissociation becomes less impor- 
tant. In the shielded regions of the cloud, the predicted 
column densities are similar to those with Go=5xl0 3 . 

The H 2 column density (or Ay) responsible of the 
warm gas observed from far-IR absorption lines (the lo- 
cation of the species in optical depth) has been tradi- 
tionally difficult to establish. A comparison between the 
observational H 2 0/OH ratio and photochemistry mod- 
els shows that the bulk of the H 2 0/OH absorption must 
arise from the surface of Sgr B2 (a maximum Ay of ~5 
to 10 mag), in agreement with the large opacities de- 
rived from the radiative transfer models. Even assuming 
an homogeneous surface cloud, most of the water vapor 
will arise from Ay<10 mag if T fc -500 K while >50% 
of water can be at Ay>10 mag if Tfc^300 K. There- 
fore, an accurate description of the thermal structure 
of the cloud (with far-IR diagnostics tracing the same 
gas) will be needed to establish more detailed conclu- 
sions. Still, in a inhomogeneous medium, several PDR- 
like clump surfaces locally heated to Tfc~300-500 K by 
low velocity shocks could be entirely responsible for the 
far-IR H 2 0/OH absorption. 

Taking into account the uncertainties implied in the 
determination of column densities in the envelope, we 
take x(H 2 O)~(l-2)xl0 -5 as a lower limit. There is a 
determination of HD column density toward Sgr B2 of 
~10 18 cm" 2 , which translates into N(H 2 ) ~10 24 cm" 2 
(Polehampton et al. 2002). However, taking into ac- 
count the huge absorption by dust at 112 /im (r ^4), it 
is unlikely that in this HD column density the cold gas, 
where high column de nsities are expected, is accounte d 
for. We derive, see l|Goicoechea fc Cernicharo 20021. 
that the H 2 column density in the warm gas is ~10 
cm" 2 . Besides, the 183.31 GHz observations imply that 
the water abundance is at least an order of magnitude 
lower in the core regions.. However, we underline the 
importance of further detection of weak H 2 lines 
and/or lines with a lesser dependence from the dust 
emission to refine the models and better constrain the 
physical parameters of the region (e.g. the temperature), 
that also determine much of the chemistry. Specially 
important will be the input of the Herschel/HIFI 
observations for water lines below ~2 THz. 

Another difficult problem is to place the origin of the 
FUV radiation field reveale d by the fine structure line 
observations in the region i|G04t) . In principle, FUV 
photons could arise from the massive stars near the 
Sgr B2(M) core and/or from another stellar population 
within the envelope itself (not resolved yet by observa- 
tions). The permeating effect of the radiation field will 
be determined by the clumpyness and inhomogeneity of 
the medium surrounding the stars, and by the energy 
of the stellar photons. In addition, X-rays observations 
could complement this scenario of widespread low veloc- 
ity shocks and UV radiation. Energetic EUV- and/or 
X-ray photons can penetrate deeper in the neutral cloud 
than FUV-photons, and thus, could also play a role 
in the chemistry, as they can also induce many photo 
ionization and photo-dissociation processes. The corre- 
lation found in the Sgr B complex between the 6.4 keV 
Fe° line and the SiO emission could also indicate that 
the X-ray sources also drive the shocks in the region 
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ijMartm-Pintado et al. 2000J) . 

Gas temperatures in XDRs can easily reach Tfc^300- 
500 K because of the more efficient gas heating by 
X-ray-induced photoelectrons from the gas (and not 
from dust grains as in PDRs). Under these condi- 
tions, neutral-neutral reactions dominate the chemistry 
and especially OH reaches large abundances. However, 
XDR models predict low H 2 O/OH<0.1 abundance ra- 
tios ijMalonev et al. 1996j) . and thus, they can not be the 
dominant scenario e xplaining the far-IR water and OH 
lines toward Sgr B2 l)Goicoechea fc Cernicharo 2002^ . 

Sgr B2 s hows diffuse emission in the Ka line of Fe° 
at 6.4 keV l|Murakami et al. 20011) . A dozen of compact 
X-ray sources have been detected within Sgr B2 cloud, 
and they may explain the whole X-ray emission found in 
the region. Many of the detected X-ray sources are not 
associated with radio (Hll regions created by massive 
stars) nor with IR sources. According to these observa- 
tions, the intrinsic X-ray luminosity toward Sgr B2(M) 
is ix^lO 35 erg s _1 , which translates into a X-ray flux 
incident on molecular gas of Fx <0. 1-0. 001 erg cm -2 s _1 
(assuming ~0.1-1 pc from the X-ray source). Hence, the 
X-ray field within Sgr B2(M) is in the low tail, or even 
weaker, than those studied by Maloney et al. (1996). 
As noted by G04, the low [Oi]63 /xm/[Cll]158 ^m 
and ([Oi]63+[Cii]158)/FIR intensity ratios observed 
in the region favor a dominant PDR or igin for these 
lines. Both shock ijDraine et al. 1983j) and XDR 
ijMalonev et al. 1996j) models predict larger ratios. 

6. SUMMARY 

We have carried out far-IR observations of sev- 
eral thermal absorption lines of water vapor toward 
Sgr B2(M) and have mapped the 183.31 GHz water line 
around the main dust condensations of the complex. The 
main conclusions of this work are the following: 

1. The detected water absorption lines are very 
opaque and arise from the warm envelope around 
Sgr B2(M). The observation of the CO J=7-6 line 
and the lack of far-IR CO lines at ISO's sensitiv- 
ities (3cr limits below ^2xl0~ 18 W cm" 2 ) imply 
that the density of such layer is n(H2)~10 cm -3 . 
The para-H20 3i3~2o2 line at 183.31 GHz shows 
~40"x40" extended emission around Sgr B2(M) 
and -40"x20" around Sgr B2(N). This is the first 
observation of that line in a GC source and repre- 
sents further evidence that water vapor is extended 
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in warm molecular clouds. 

2. LVG and nonlocal radiative transfer calculations 
have been carried out to extract the water va- 
por abundance and to constrain the physical pa- 
rameters of the absorbing/emitting regions. Be- 
cause of the radiative excitation by dust photons, 
the far-IR water lines are not very sensitive to 
Tfe. Taking into account the analysis of the re - 
lated species OH l)Goicoechea fc Cernicharo~2 002'l. 
the water absorption must arise from warm gas 
at similar temperatures, i.e., from 300 to 500 K. 
For this warm envelope, we found A r (H20)=(9 ± 
3) x 10 16 cm~ 2 . An important fraction of the 
183.31 GHz emission arises from the inner, denser 
and colder gas located closer to the main cores. We 
estimate a water abundance of a fewxlO -7 in the 
denser regions. 

3. Photochemistry models show that a component of 
warm gas, ^300-500 K, is needed to activate the 
neutral-neutral reactions and reproduce the large 
H 2 and OH column densities observed in the en- 
velope. We show that OH and H2O far-IR obser- 
vations toward Sgr B2 are surface tracers of the 
cloud (a maximum Ay of 5 to 10 mag). We found 
X(H 2 0)~(l-2) x 10~ 5 in these regions. Although ir- 
radiated by FUV, and possibly more energetic pho- 
tons, affecting the H 2 0/OH ratio in the outermost 
layers, a clumpy structure for the PDR is needed. 
Alternatively, low velocity shocks could maintain 
the gas heating through the envelope. 

Due to the complexity of Sgr B2 (and also of the GC 
ISM as a whole) a multiple scenario is needed to ex- 
plain the modest angular resolution far-IR observations. 
The input of chemical models, and higher sensitivity and 
larger spatial resolution far-IR observations will lead to 
a better understanding of the GC environment. 
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TABLE 1 

Summary of H2O spectroscopical and observational data presented in this work. Uncertainities refer to 

FITTING ERRORS. THE VALUES PROVIDED FOR THE 183 GHz LINE CORRESPOND TO THE THE DIRECTION OF SGR B2(M). 
IN THIS CASE, THE AREA OF THE OBSERVED LINE IN K KM S _1 IS PROVIDED INSTEAD OF THE FLUX. NOTE THAT ENERGY 
LEVELS TAKE INTO ACCOUNT THAT BOTH ortho AND para GROUND-STATES ARE AT K. 



species 


transition 


A(/im) 


i/(GHz) 


Eupper (K) 


A / — 1 \ 


F; (W cm -2 ) 




F/F c 


Av (km s _1 ) 


P-H2O 


3l3 


— 220 


1635.439 


183.3 


205 


3.53e— 06 


1120±60 K km 


s- 1 




30±3 


o-HfO 


2l2 


- I01 


181.053 


1655.8 


79 


5.45e-02 


-(2.68±0.26)e- 


-18 


0.73±0.02 


52±2 


0-H2O 


221 


-2l2 


180.488 


1661.0 


160 


2.99e-02 


-(3.44±0.46)e- 


-18 


0.74±0.01 


64±3 


0-H2O 


2l2 


- I01 


179.527 


1669.9 


80 


5.47e-02 


-(6.76±0.29)e- 


-14 


0.06±0.03 


204±23 


0-H2O 


303 


-2l2 


174.626 


1716.7 


163 


4.94e-02 


-(1.86±0.46)e- 


-18 


0.82±0.02 


44±2 


P-H2O 


3l3 


— 2()2 


138.527 


2164.1 


205 


1.22e-01 


-(1.27±0.03)e- 


-18 


0.89±0.02 


41±2 


0-H2O 


432 


- 4 2 3 


121.719 


2462.9 


516 


1.20e-01 


-(5.40±0.29)c- 


-18 


0.96±0.01 


62±9 


0-H2O 


221 


- I10 


108.073 


2773.9 


160 


2.52e-01 


-(6.07±0.05)e- 


-18 


0.64±0.02 


53±2 


p-H^O 


220 


- Ill 


102.008 


2938.9 


194 


2.52e-01 


-(9.38±2.62)e- 


-19 


0.96±0.01 


72±8 


P-H2O 


220 


- Ill 


100.983 


2968.7 


196 


2.55c-01 


-(6.02±0.11)e- 


-18 


0.66±0.01 


55±9 


P-H2O 


322 


-2n 


89.988 


3331.4 


297 


3.45c-01 


-(1.99±0.14)e- 


-18 


0.90±0.01 


67±2 


0-H2O 


321 


-2l2 


75.380 


3977.0 


271 


3.25c-01 


-(4.32±0.07)e- 


-18 


0.69±0.02 


51±2 


P-H2O 


331 


— 220 


67.089 


4468.6 


410 


1.20e+00 


-(1.36±0.08)e- 


-18 


0.93±0.03 


76±8 


0-H2O 


330 


-2 2 i 


66.437 


4512.4 


376 


1.22e+00 


-(8.39±0.14)e- 


-19 


0.92±0.02 


46±3 



TABLE 2 

Summary of far-IR CO spectroscopical and observational data presented in this 

WORK. NO FAR-IR CO LINE HAS BEEN DETECTED. 3cr UPPER LIMIT LINE FLUXES ARE 

TABULATED. 



species 


transition 


A(^im) 


iy(GHz) 


Eupper 


Aij (s 


- 1 ) 


F ; (W cm" 2 ) 


F/F c 


CO 


14- 13 


185.9 


1611.8 


581 


2.95c- 


-04 


<2.28e-18 


<1.007 


CO 


15 - 14 


173.6 


1726.6 


663 


3.64e- 


-04 


<1.02c-18 


<1.003 


CO 


16 - 15 


162.8 


1841.3 


752 


4.42e- 


-04 


<1.14e-18 


<1.003 


CO 


17- 16 


153.3 


1956.0 


846 


5.31e- 


-04 


<1.49e-18 


<1.003 



